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THEORY OF SPECTRAL RADIANCE OP POLLUTANTS AT SEA 


This report concerns remote measurement of pollutants dumped 
in the sea, not oil slicks, but soluble pollutants that change 
the color of the water. The sensor is a spectral radiometer 
that flies over the polluted area and compares its spectral 
radiance (color) to that of surrounding clean seawater. It 
is fairly easy to make a qualitative determination that a 
plume of pollutant occupies a certain area. However, the 
goal liere is quantitative, namely to infer the concentration 
of pollutants using the measured radiance of the sea compared 
to laboratory measurements of reflection and transmission 
spectra of the pollutants. Many previous papers and reports 
have discussed remote sensing of ocean color. The main 
contribution here is how to uses data quantitatively that 
compare the polluted area to surrounding clean water. 

We treat the subject in three steps. The first section defines 
the quantities involved and describes means for measuring 
them. The second section derives the equations for remote 
sensing with a low-flying aircraft, in which case the absorption 
and radiance of intervening air is negligible. The third 
section applies to high-flying aircraft and satellites, in 
which case the radiance of intervening air is the major problem. 
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1.0 


FUNDAMENTALS 


This section describes measurements of light intensity, reflection, 
and transmission. 

1.1 Irradiance 

The brightness of the sea surface is proportional to the amount 

of incident sunlight measured as an irradiance, so much power 

2 

(watts) per unit area (m ) . This incident light is a combination 
of collimated direct sunlight and diffuse sky light. The angles 
of incidence are not too important as long as the sun is fairly 
high in the sky because refraction at the sea surface bends the 
rays toward the nadir. Fig. 1, so that all the light is directed 
more or less downward, 

A meter that measures irradiance (irradiometer ) is nothing more 
than a lensless light detector with a precise area A of detection 
defined by a mask, Pig. 2. The detector measures power P, and 
irradiance is given by 

H = P/A (1) 

Since A is constant, a practical meter is calibrated directly in 

2 

irradiance units, such as watts/m . A suitable antireflection 
coating on the detector helps to make the meter insensitive to 
the angle of incidence. To measure light incident on the sea, 
it is important that the irradiometer be exactly horizontal. 
Obviously it will intercept more sunlight if tilted toward the 
sun. 

Since we are interested in color of the sea, we need a spectral 

. 2 

irradiometer which measures spectral irradiance, watts/m /nm. 

The irradiometer is converted to a spectral irradiometer by 
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covering it with a series of absorbing filters of bandwidth B. 
(Interference filters would be sensitive to angle.) Then the 
spectral irradiance is given by 

= H/B = P/AB (2) 

Numerous measurements of irradiance have been made and compiled 
in texts on ocean optics, e.g. Jerlov [Ref. 1] pp 71, 128 (Figs. 
63 and 65b) , and 144 and Jerlov and Nielsen [Ref. 2] pp 88 and 
39 9. Two sunmiaries of irradiance appear in Figs. 3 and 4. Note 
that just beneath the surface, downwelling irradiance in the 
blue-green is about 20 times the upwelling. This ratio increases 
to 35 in the red. In a later estimate we use the value 

H(510 nm) =1.2 W/m^/nm 
1 . 2 Radiance 

The remote sensor that looks at the sea cannot be a lensless 
irradiometer . An imaging sensor is required to distinguish 
a plume of polluted water from surrounding clean water. The 
amount of light power this sensor detects from each pixel is 
proportional to the solid angle of rays that it intercepts, 
and so the appropriate measure of sea brightness is power 
(watts) per unit area (m**) per unit solid angle (steradian) , 
a quantity called radiance, L. 

Figure 5 shows the optics in a radiometer. Power received at 
the detector is given by 

_ /area of resolution cell\ /solid angle intercepted 

P ~ L X ^Fig. 5, lower right / \Fig. 5, lower left 

From the geometry, these quantities are evaluated as follows; ’ 
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m nm 



Fig. 3. Comparison between spectral distribution 

of downward ( — ) and upward ( ) irradiance 

for a solar elevation of 55-60* in the 
Sargasso Sea. (Lundgren and Hojerslev, 
Kobenhavns Universitet, Rep. Inst. Fys. 
Oceanog. 14 (1971)) 



/t£LAri^£ £N£RGY 



n 

3 

O 

OQ 


a 

o 

c 

3 

€ 

ra 


3 

>Q 


H* 

►1 

01 

& 

H- 

0) 

3 

O 

o 


s: 

\ 

3 

K> 

\ 

3 

3 


Fig. 4. Compilation of irradiance data by G.L. Clarke 
and G.C. Ewing, Chapter 17 in Ref. 2. Spectra 
of the energy of upwelling light on 13 August 
1969 in Buzzards Bay at the altitudes shown 
(measured by TRW spectrometer) and at 0*20 m 
below the sea surface (measured by Scripps 
spectrometer) as compared with downwelling 
irradiance of the NASA Standard spectrxim 
outside the atmosphere, and the downwelling 
irradiance measured on 14 August 1969 at the 
sea surface. For the "above surface" curve 
see text. The energy scale on the right is 
for downwelling light; that on the left is 
in relative units, except that the values on 
the left also represent lO'^w/m^/nm for the 
upwelling light at 0*20 m. 
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P - L 



Aa 
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where A - ttd 2/4 is the collecting aperture , and D its diameter, 

f the focal length, and h the altitude of the sensor. Note 

that h drops out, and Aa/f is a constant of the optics. This 

means that the detector output can be calibrated directly in 

2 

terms of radiance, watts/m /nm. Also note that the radiance 
of the image is the same as that of the object, neglecting 
absorption in the lens; 


*‘i^ 


(A/f^)a 


111 

Aa 


(4) 


Again we are interested in color of the sea, a.nd so spectral 
radiance is required: 

L;^ = L/B (5) 

Absorbing filters could be used, but in this case the optics 
can be used to collimate light in a suitable way tc use a 
grating or prism monochromator. These details are beyond the 
scope of this report. From this point on, we drop the subscript 
X with the understanding that all radiances and irradiances 
are spectral. 


Austin [Ref. 3] has compiled graphs of radiance that are 
particularly useful for remote sensing. Four of these appear 
in Fig. 6 and 8 through 10. Fig. 7 defines the notation ha 
uses. The following features in these date influence the 
techniques discussed later: 

► L (510 nm) = 0.9 x 10 w/m /sr/nm regardless of 
the color of the water (Fig. 6) 

^ In the deep red, sky radiance reflected by the 

surface (L ) exceeds that which has penetrated 

the water (L ) . The crossover falls in the range 
w 

580 - 650 nm. 
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Spectral radiance, W/m /sr/nm 




IVAVEIENGIH (nm) 

Fig, 6. Upwelling spectral radiance just below 
the surface computed from measurements 
at depths of 2.2 and 7.7 m. 
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^ From high altitude or outer space, the observed 

radiance (L ) is dominated by the veiling luminance 
of the intervening air (L*) . 


1.3 Diffuse Reflectance 


Consider direct sunlight incident on a rough surface as shown 
in Fig. 11. The ratio of reflected to incident power is 
reflectance p, a dimensionless quantity. In order to relate 
L seen by the radiometer to H, the incident light, we need 
the ratio L to H, a quantity with dimensions of per steradian : 


L 

H 



( 6 ) 


where is 

is scattered. 


the effective solid angle over which the light 
We call this ratio the diffuse reflectance.* 


Consider the bulk reflectance R of a liquid, excluding reflection 
at the surface, which we treat separately. This liquid can 


* In the case of a perfectly white diffuse reflector, p = 1 
and "n"/ the effective solid angle for a lambertian 

(cosine law) surface. Hence = 1/n. This relationship 

provides a means for calibrating a radiometer against an 
irradiorneter given a near-perfect surface. For airborne 
instruments, a cloud bank serves this purpose very well 
providing it is thick enough to stop direct sunlight (cannot 
see the sun from below) but not thick enough to be an absorber 
(dark thunderhead) . In this case 

V "-down “ . <’> 

where values of L are measured from above the clouds looking 

down, and below looking up and for, a cloud of the right thickness 

they are about, equal. 
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be either clean or polluted seawater; in either case it contains 
particulate matter that reflects light diffusely. It is 
possible in principle to derive an expression for R using the 
attenuation o of the medium and the diffuse reflectance adz, 

Fig. 12, caused by particulate matter {plus Rayleigh scattering) 
in each infinitesimal layer. A simple expi'easion useful later 
applies to tlie case in which the absorption rate, a, dominates 
over scattering; 

a > s = a(0)duj or a > 4ttc;(180°) 

Multiple scattering is negligible in this case, and the single 
scatter i-ays in Fig. 12 lead to the expression 

„ /•«> -az , j . -az /*“„“2ciZj 

^ - Jo ® (adz)e = a e dz 

R = a/2a, a > 4Tfa(180'’) (8) 

It is possible to treat more general cases, but we see no point 
in relating R to the bulk optical properties of the liquid. 
Instead it is more reliable simply to measure reflectance in 
the laboratory for various concentrations of the polluted 
layer. This is especially true in view of the fact that some 
pollutants undergo chemical change as the concentration changes . 
For example, the acid waste that DuPont dumps precipitates iron 
when diluted. 

1.3.1 Reflectance of Pollutants 

Figure 13 shows an apparatus suitable for measuring diffuse 
reflectance. A known irradiance H illuminates a cylindrical 
column of liquid (thickness z) . .A radiometer in a side arm 
measures L, and then 

RT = L/H (9) 
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radiometer 






Fig. 13. Apparatus for measuring 
diffuse reflectance of 
pollutants 


Fig. 12. Single scatter 
event 



Here T is a transmission factor for the air/liquid dielectric 
interface (discussed later) . It is just a well to measure the 
product RT instead of each separately, because the product is 
what appears later in the equations for remote sensing. 

Several features of this apparatus (Pig. 13) are noteworthy. 
First the cylindrical wall below the surface should be highly 
reflecting, silvered preferably on the inside (and coated to 
prevent corrosion) . This simulates the effect of an unbounded 
ocean by replacing light that would come from outside the 
cylinder (dashed ray lower right) by; a 'ref lected ray on the 
inside. Other surfaces should be black to reduce stray light. 
To avoid edge effects, the incident beam must overfill the 
cylinder, and the radiometer's field of view must not extend 
to the side walls. 


1.3.2 Diffuse Spectral Reflectance (Color) of Clean Seawater 


We tend to think that the color of seawater is obvious , s '^mewhere 
between blue and green, but actually the apparent color is 
contaminated with a lot of sky reflection, somewhere between blue 
and white, and the correction is not always straightforward. 
Probably the only reliable means to measure the true color of 
seawater alone, i.e. R^, its spectrum of diffuse reflectance, 
in all weather requires a boat trip to the site in question. 

One measures radiance (upwelling) with a radiometer just 
under the surface looking down, and irradiance (downwelling) 
with an irradiometer just under the surface looking up. Then 


\ = 


In a few sites where pollutants are dumped regularly, it may be 
practical to do this enough times to establish the variations 
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of R with season, hour of the day, and sea state. Then, 
from that time on, R^ may be considered a known spectrum 
in the remote sensing equation. 

Obviously we need a fall-back position for areas that do not 
justify all that labor. One alternative is to fly over the 
area.s and take measurements when the sea is calm and the sky 
uniform. Under these circumstances, the remote sensing 
equation (Eq. 17, discussed later) may be solved for R^^ since 
all other quantities are either known or easily measurable. 

One has to hope that the values of R^ found in this way also 
hold when the weather changes . 

The final fall-back position is to use only two colors where 
R^ is fairly invariant from place to place and therefore 
known. These values are 

Green: R^(510 nm) = 8 10 ^ sr ^ 

(1 

Near infrared: R^(700 to 1100 nm) = 0 

These numbers require some explanation. The green color is 

where blue water, blue-green water, and green water all have 

the same reflectance. Fig. 6. This assumes that the water in 

question is moderately deep and its color is not influenced 

by bottom reflection, silt, or gelbstoff. Unfortunately we 

know of no case in which R (510 nm) or any other wavelength 

w 

has been measured directly. We calculated the value above 
from a ratio L/H in which L and H were typical values measured 
in widely different parts of the world by different groups 
who presumably do not compare instruments. We took L from 
Pig. 6 and H is an average from Figs. 3, 4, and others like 



them. This situation should be corrected because any quantity 
that is fairly invariant should be pinned down and published 
as a universal semi-constant whose range of variation is known. 

The second of Eqs . 10, R^^(700 to 1100) = 0, simply means 
that R is negligible. To estimate the actual value, let us 
use Eq. S with 

-3 

a(160°) = 10 per meter per steradran 
a = a = 0.5 per meter. 

The former avoids vT fairly sharp peak at ISO® and represents 
a fairly large value of the scattering function in a backward 
direction. The value of the attenuation constant a comes 
from Fig. 14. Now Eq.8 gives 

_3 

R = a/2a = 10 per steradian, 

an order of magnitude less than the green value, which already 
represents a very dark scattering medium. (Deep seawater 
looks almost black compared to most diffuse reflectors.) The 
useful range of infrared cuts off at about 1100 nm where light 
penetrates only about a centimeter of water. Fig, 15, hardly 
enough to see an ’appreciable amount of diluted pollutant. 

In view of the difficulties in determining R^, it is advisable 
to formulate the remote sensing equations in such a way that 
this quantity drops out. This is done in Sec. 2, but the 
result is not usable in all circumstances. 

1.4 Sky Reflection 

Reflections from fog, foam, spray and particulates in the sea • 
are all of the diffuse type described in the preceding section. 
There is one exception, reflection at a calm sea surface. 
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This is a simple specular reflection with no angular spread/ 
and so collimated direct sunlight is still collimated after 
reflection, and the radiometer sees a reflected image of the 
sun. The brilliance of this image is overpowering and masks 
the color of the sea. Any radiometer looking down at the sea 
must be pointed away from the sun's image to obtain useful 
data. With this precaution, the downlooking radiometer sees 
the color of the sea, but it is still contaminated by a 
significant amount of reflected sky radiance: 


= rS 


where S is the radiance of the sky as measured looking up at 

the angle supplementary to the direction of as Fig. 16 

shows. The extent of this contamination is considerable: 

L is about 1/3 of the desired radiance L from the water 
r w 

(Figs. 8, 9, and 10) throughout most of the spectrum, but 

the two cross over in red and the contamination exceeds the 

signal. Since no angular spreading occurs at reflection, 

2 

and both and S have spectral radiance units (W/m /sr/nm) , 
it follows that r is a dimensionless number like the reflectivity 
of any mirror surface. 


Values of r are plotted in Fig. 17, which shows that it is 
dependent on polarization, the symbol "j." referring to polarization 
(E-vector) perpendicular to the plane of incidence and "ii" to 
polarization in the plane of incidence. Since skylight is 
also polarized, the equation above becomes two equations: 









( 11 ) 


For small angles, i.e. directions within about 20° of nadir. 
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Fig. 19. Sky luminance distributions 
in sun-zenith plane for 
several locations and solar 
zenith angles, 65 . Data 
compiled by Austin [Ref. 3] 
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the reflectivity is given by 



the well-known Fresnel formula for dielectric reflection at 

normal incidence in which n = 1.341 is the refractive index 

w 

of seawater and n„ = 1.0 is that of air. 

a 

At other angles r is given by formulas that are straightforward 
but complicated and are given in numerous physics and optics 
texts. They are of little if any practical value here because 
the sea surface is ruffled by wind which averages the reflectivity 
over a range of angles. Gordon [Ref. 4] has computed the average 
for various wind speeds with the results given in Pig. 18. 

The wind-ruffled sea surface blurs the reflected image of the 
sun into a glittering pattern called sun glint. The radiometer 
should be pointed away from the sun to avoid this glint as 
much as possible, Fig. 16. However, to reduce sky contamination 
to a minimum, the radiometer should look somewhat farther from 
the reflected sun angle, as Fig. 19 shows, about 75° when the 
sun is high in the sky to 90° when it is low. 


By coincidence, this same direction gives the maximum polarization 
of reflected skylight in the perpendicular direction due to both 
factors in Eq. 11: 


> r and S, > S 
ti i 


II 


The former fact is apparent from Fig. 17 and the latter results 
from the dipole radiation pattern of molecular scatter shown 
in Pig. 16. Unfortunately Fig. 18 is computed for unpolarized' 
light and we do not have a similar 'one for the separate components 



When the sky is mostly clear, so that light has the directional 
qualities shown in Fig. 16, then the radiometer should be f'ointed 
to avoid reflected skylight and a polarization analyzer should be 
used to rejecb the perpendicular component. I'he radiometer's 
field of view should be centered about a direction 75 to 90® 
away from the reflected sun image in the plane containing that 
image and the nadir. This technique gives four advantages: 

^ Avoids sun glint 

^ Rejects sky contamination by looking at the darkest 
part 

|> Further rejects skylight by taking advantage of its 
polarization 

Takes further advantage of the polarization by 
Fresnel reflection. 

If the sky is clear (so that molecular scattering dominates 
over particle scattering) , then a polarization filter almost 
completely eliminates sky reflection at the longer wavelengths.* 
This is fortunate, because it normally causes most of the 
veiling luminance at the longer wavelengths; compare and 
in Figs . 8 , 9 , and 10 . 


*Some blue and violet remain, which have been depolarized by 

double scattering in the atmosphere, but this falls off as 
-8 -4 

X , a factor of X for each of two Rayleigh scatters. This 
is why the sky looks violet in some photos taken through 
Polaroid at right angles to the sun. 
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1.5 


Transmission 


We shall use transmission factors T to account for loss of 
light in propagating through lossy media; 

T - atmosphere {satellite case) 

T “ sea surface 
s 

Tp - polluted layer 

In every case, light that passes down through a layer must 
return through the same layer to be detected, and so we use 
capital T to represent round-trip propagation. However, this 
does not imply that the loss is the same each way. For example, 
when the sun is rather far off the zenith. Fig. 20, and the 
radiometer is looking straight down, the transmittance through 

some layer going down is T = e ' ^ ' while going up it 

is considerably less, namely T = hence, 

T = +?■ = exp [-aza (1 + sec 6)] (13) 

t 

In estimating T, light that is scattered through small angles 

should not be counted i.-s lost. A cutoff angle is not easy to 

define, and will have to be determined by empirical matches 
to the remote sensing equations using data taken at sea. Fig. 21 
shows a summary of estimates of the unscattered transmittance of 
the atmosphere needed for the case of a. satellite sensor. 

1.5.1 Polluting Layer 

For measuring Vt" a suitable absorption cell is shown in Fig. 22. 
The inside cylindrical wall is silvered for the same reason 
discussed in connection with Fig. 13. Similarly, the known 
incident irradiance overfills the stop aperture, and the 
detector is large enough to catch all light at the output. 
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radiometer 



Fig.' 7 . 0 » Round-trip attenuation 



Fig. 21. Vertical beam transmittance of the atmospheric 
models (curves A, B, C, D) and as measured 
by Guttman (curves a and b) . Applied Optics 7 
(1968) p. 2377 
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Transmission varies approximately exponentially with distance z, 




but not exactly owing to multiple scattering in thick layers. 
The attenuation rate may be expressed 


a = a + a C 
w o 

where a is the attenuation in the water that dilutes the 
w 

pollutant, ct^ is attenuation of undiluted pollutant, and 
C is its concentration. Con±>ining these equations gives 


The factor Cz in parentheses is the area concentration of 
pollutant, so many grams per square meter, which says in effect 
that it makes no difference how the pollutant is distributed 
in depth. This is true so long as. there is no chemical change 
that depends on concentration, as there is in 'the case of the 
DuPont acid dump. In the latter case, is a function of C. 

It is prudent not to assume that Eq. 14 holds over any wide 
range of C or z until verified in the laboratory. Over the 
intervals where it does hold, Eq. 14 provides an interpolation 
formula. 


1.5.2 


Surface Transmission 


Near vertical transmission through the- sea surface is given by 

2 / „ \ 2 


( 4n n \'^ / n \ 

- \ / ^ 1 

<"a "w> / \ «/ 


(15) 


The first factor is round-trap Fresnel reflection at the dielectric 

2 

interface ((1-r) , see Eq. 12). The second factor represents the 
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reduction of radiance rising out of the water due to angular 
spreading, Fig. 1. Using n^ = 1.0, n^ = 1.34 as before gives 

Tg = (0.96) (0.56) = 53% (16) 

The fact that the first factor has little effect implies that 
surface roughness is not a big factor in determining , as 
shown in Fig. 18, except, that patches of foam cause significant 
changes by scattering light upward and sideways. 
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2.0 


REMOTE SENSING EQUATIONS 


This part of the report derives equations for determining 
pollutant concentration in two cases. The first is for a 
low altitude aircraft, and the second for high altitude, most 
likely a satellite. 

2 . 1 Low Altitude 

The significance of this case is that low altitude eliminates 
all atmosphere effects except reflected sky radiance, and we 
can deal with that by quite effective means. As discussed in 
Sec. 1.4, a good choice of aspect angle plus a polarization 
filter eliminates most of the interference from the sky when 
it is clear. At other times it is necessary to measure the 
spectral radiance of the sky directly with a sky radiometer 
so that the reflected radiance may be subtracted out of the 
sea data, Fig.. 23. The best arrangement is to combine sea 
and sky radiometers into one instrument with the same detector 
and amplifier, so that any fluctuations in sensitivity affect 
both measurements equally. An optical switching or scanning 
arrangement causes the detector to look at the sky pericdically 
between looks at the sea. (It need not look at the sky very 
often.) In any case, the pointing of the sky and sea radiometers 
must be at supplementary angles to the nadir as indicated by 
the reflection geometry in Pig. 23. . 

When the sea has normal roughness, the reflected image of the 
sky is blurred. This causes no problem when the sky is uniform, 
either clear or overcast. However, when the sky has broken 
clouds, the optics in the sky radiometer must blur the direct 
image in a manner that simulates the blurred reflected image 
that the sea radiometer sees. The proper amount of blur may 
be determined by measuring the angular diameters of the eliptical 
swn glint pattern. 
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Quantities fo 
sensing equat 
altitude case 


Polluted layer 
over clean 
seawater 




The radiance seen over clean seawater is given by 





^c = "^^s\ 


( 17 ) 


As Fig. 23 shows, the first term represents light that penetrates 
the water, and emerges, and the second term light that reflects 
from the surface. The latter is L^, Eq. 11, either S,, or 
r,iS|, + depending on whether a polarization filter is 

used. The first term in Eq. 17 is the larger by a factor 
of about three through most of the visible spectrum, but the 
two cross over near the red end, Figs. 8 through 10. Recall 
that every quantity in Eq. 17 is a spectrum. Their colors 
are as follows: 

- blue’ to green (water color) 


Also recall that Eq. 17 implies three spectral measurements of 

light intensity. The radiance of the sea and sky, L and S 

o 

just discussed, and the irradiance H. 


j blue when clear 
(white when overcast 


H 

T. 


R 


white 


w 


S - 


Now consider the changes that occur when the water is polluted. 
An oil slick changes . r considerably, but this discussion is 
restricted to the case of soluble pollutants, which have 
negligible effect on r and merely change the color of the 
water. Fig. 24 shows, the reflection and transmission factors 
that lead to the following equation:- 


L = HT (R + T R ) + rS 
P s ' p p w' 


(18) 


The radiometer flying over a plume of pollution measures both 
L and L while also monitoring H and S. With these results. 


Eqs . 17 and 18 are two equations to be solved for two unknowns. 

One unknown is the desired pollutant concentration implicit 

in its reflectance R and transmittance T . The other unknown 

P P 

can be one of the quantities R (water color) , T (surface 

5 

transmission) , or rS (reflected sky) , that is not easy to 
measure or estimate accurately. There are different cases 
in which any one of these may be poorly known. However, we 
postpone that discussion to complete the mathematics of the 
remote sensing equations. 

When the color of the water R is unknown and eliminated 

w 

from Eqs. 17 and 18, the result is Eq. 19 in the table on 
the following page. When surface transmission T is unknown, 
the result is Eq. 20. Tnis equation has the added advantage 
that irradiance H is eliminated and therefore does not have 
to be measured. When reflected sky rS is unknown, the result 
is Eq. 21. In general one must solve the appropriate remote . 
sensing equation 19, 20, or 21, by -an iterative algorithm: 
assume a pollutant concentration, test the corresponding 
R (A) and T (A) in the equation, and evaluate the fit. Assume 

•i •> 

another concentration and evaluate, and so on to find a best fit. 

A limiting case of interest occurs when the pollutant is a 

clear ‘ colored liquid, i.e. R_ = 0. In this case the remote 

P 

sensing equations (19 - 21) become Eqs. 22 and 23 in the table. 

Whenever L > L , this case cannot apply because the transmittance 
pc , ■ 

T has to be less than one. Eg. 22 is particularly nice because 
P 

so many quantities drop out and therefore need not be measured. 

Another limiting case is T^ = 1, i.e. the pollutant does not 
appreciably shade the water beneath it.* This gives Eqs. 24 and 25 

*Strictly speaking, if T =1, then R = 0, but reflectances 

P P 

at sea are very small, so R can be .3% per steradian and T = 95%. 

P p 

is close enough to 1 . 
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ABLE OF REMOTE SENSING EQUATIONS 
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Equation 24 is a nice case because so few quantities are 
needed . 

Let us return to the postponed questioi., which quantity should 
be treated as the second unknown? Much experimental work is 
needed to establish the circumstances reliably, but some 
general guidelines are as follows : 

When the sea is fairly calm treat R,. (water color) as unknown 

and use Eq. 19 . Sec. 1.3.2 discussed the fact that the color 

of seav/ater is not as obvious as it seems, and in fact measurements 

in situ are required to reliably determine Accordingly, 

R should be treated as the unknown whenever the sea is calm 
w 

enough to be free of whitecaps so that and r are fairly 
well known. 

W hen the sea is choppy with foam and spray and the sky is 

uniform, assume T „ is unknown and use Eq. 20 . One must work 

in a location where R^ has previously been established or else 

use only two colors, green and infrared, as discussed in 

Sec. 1.3.2. Foam and spray make it impossible to estimate 

either T or r. However, r may be measured indirectly leaving 
s 

T as the unknown. For X ~ 700 nm, T is negligible, so Eq. 17 
s s 

reduces to 

L^(X > 700) = rS(X > 700) 

The sky radiometer measured S directly-, and sor{X > 700) = 

L /S. Moreover, the dependence of r on X is very weak and 
c 

known, so r is established throughout the spectrum leaving 

T as the unknown, 
s 


35 


When the sky and sea are too complex to find a reasonable 


estimate of blurred reflected sky radiance, treat rS as the 


unknown and use Sc 


When both T and rS are indeterminate, 
s 


rS should be treated as the unknown in the equations, because 
Tg is never too far from the factor (n^/n^) for the spread 
in solid angle. Other losses at the surface are not large 
(but not negligible either) . 

2,2 High-Altitude Sensor, Especially Satellite 


Let and denote the spectral radiance (clean and polluted) 
measured by a sensor at high altitude. These quantities are 
related to the low-altitude measurements discussed in Sec. 2.1 
as follows: 

i = L T + L* 1 

c c a a . 

L = L T *i* L* I 
p p a a I 

Here is the one-way upward transparency of the air path 

a 

and L* the' radiance (blue) of the air along the path. In 
other words, the luminance w'e wish to observe is attenuated 
by T and veiled by L*. Unfortunately the path radiance L* 
is about seven times larger than the terms with the desired ‘ 
information. This means that L* must be subtracted off with 
considerable precision to avoid contaminating the ocean data 
with its errors. 


Recall that and in these equations are given by Eqs 
17 and 18: 


L = HT R + rS 
c s w 

L_ HT_(R_ + T_R.) + rS 




In Sec. 2.1 we processed these two equations in three different 
ways appropriate to three different cases. However, in the 
high-altitude case, we always take the difference: 





c s a 




- 


(27) 


This not only eliminates the largest term L* (and its uncertainty ) , 
but also rs in which neither r nor S is known very well from 
high altitude. (Sea state and low-lying hydrometeors may be 
in doubt.) 


The factor H in Eq. 27 is irradiance at the sea surface, a 

quantity not normally available when the only measuring platform 

is at high altitude. Therefore, we replace it by the solar 

irradiance H measured at high altitude times some estimate 
s 

of downward atmospheric transmission; 


H = T, 
■ . s a 


^L = (T Ta) T^ 
s a s 


(R + T R - R ) 
p p w w 


It may seem that the product +a t should be replaced by T 

3 , a 

to represent round-trip propagation as in Eq. 13 and other 
places. However, there is a subtle distinction between the 
two directions that needs to be emphasised. In the downward 
direction scattered light counts as long as it reaches the 
surface and thereby contributes to illumination. In the 
upward direction, this is not the case. As in Fig. 25, the 
polluted area is usually quite small compared with the altitude 
of the sensor. As a result the scattered light is diluted 
and completely lost in the scattered light from surrounding 
clean water. A distinction between polluted and clean areas 
is possible only in the unscattered light because the radiometer 
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images them and compares radiance on each side of the image 
boundary. Thus light scattered in upward propagation drops 
out of the difference (Eq. 27) but unscattered light contributes, 
and we emphasize this with a chage of notation: 

D ~ T including scattered light (D stands for down 
and diffuse) 

U - excluding scattered light (U stands for up 
and unscattered) 

Finally we have the remote sensing equation; 

(Rp - (1 - Tp) R„] (28) 

The evaluation of this expression is not nearly as satisfactory 

as the corresponding equation for low altitudes (choicsi of 

Eq. 19, 20, or 21) . The direct measurements of H (at the 

surface) and S (sky) in the low altitude case have been 

replaced by quantites that involve indirect evidence and even 

conjecture. Clouds near the line of sight affect Foam, 

spray, and possibly other low-lying hydroraeteors affect T . 

s 

One can estimate a correction based on whiteness of the reflection 
from clear areas. However, it cannot be very accurate because 
one does not know how much of the whiteness to attribute to 
clouds (reduce D^) , surface foam and spray (reduces T^) , ,or 
to reflected sky (clouds), which has dropped out of Eq. 28. 
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